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A technique for avoiding infinite integrals in the calculation of the one-loop diagram con- 
tribution to the vacuum polarization component of an atomic energy level is presented. 
This makes renormalization unnecessary. Infinite integrals do not occur because, as it is 
shown, no delta functions are required for the Green's functions. Thus there are none to 
overlap. This procedure is shown to produce the same formula as the one obtained by 
^ dimensional renormalization. 
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X) I. INTRODUCTION 

D . 

^ ' Divergent integrals appear in the calculations of the Lamb shift in a number of papers 

^ . and text books. They are removed by one of the processes called "renormalization" . The 
I present author once asked Richard Feynman if these infinite integrals represent a defect 

■ in the fundamental equations of QED or do they originate in the method used to solve 
^ I the equations. He paused for a moment and then replied, "I don't know". The same 

question to Paul Dirac produced the answer, "The basic equations are at fault." We shall 
Q i see evidence here, however, that there may be no defect in the fundamental equations, but 

■ that the infinities may arise from the use of faulty Green's functions as they do in many 
5 ! one-loop calculations. 

O ■ Actually, this is a solution to a very old problem that held up the development of 

^ . quantum electrodynamics (QED) in the years preceding 1947. Feynman^ noted that the 
I I infinite integrals that were blocking the calculations in those days were present because of 
the overlapping of the delta functions in the photon or fermion Green's functions. He and 
a number of others have invented forms of renormalization to subtract the infinite terms. 
However, they were unable to avoid the appearance of these infinities in the first place. A 
• ^ , means for avoiding the infinity in the electron self-energy has been published"^. However, 
rS I the renormalization procedure for the vacuum polarization contribution to the Lamb shift 
. absorbs an infinite term into the electric charge. No way of guaranteeing a finite charge 
" " ' from this procedure has been suggested. A way of avoiding this problem will be shown 
here. 

We note that Green's functions are often used to solve problems in potential theory 



involving Poisson's equation. These Green's functions often have the form, |r2 — 
This is to be compared with Feynman's"^ propagator which is. 
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K+(2,l)=z(z^ + m)/+(2,l) (1) 

where 

/+(2, 1) = -{4n)-^S{s'^) + {m/Ms)Hf\ms) (2) 

where s = (t^ _ x2)V2 for t^ > and s = -i{x^ - t2)l/2 for f < x^. ifj^^ is the 
Hankel function. The three dimensional Green's function used in potential theory has no 
delta function while the four dimensional Green's functions for the electron and photon 
do. When we consider the vacuum polarization electron-positron loops, it will be clear 



that the delta functions from the two fermion propagators wiU have the same argument. 
This is an example of a case where an infinity occurs. 

In Section II, a new fermion Green's function will be derived in a manner similar to the 
derivations often used for potential theory. In Section III, they will be applied to obtain 
the formula for the contribution of the fundamental single-loop vacuum polarization to the 
Lamb shift. This formula is identical to that derived using dimensional renormalization. 
Some implications of this result will be discussed in Section IV. 

II. THE FERMION GREEN'S FUNCTION 



In order to ease the visualization of the steps required to define a covariant Green's 
function, we make use of a Wick rotation to obtain a complex time coordinate defined by 
a;^ = ict and a 7^ = i-f^ so that p = 7^94 + ^^dj. We define to equal {x^f + (a;2)2 + 
(x^f + (x^) and = (^4)2 + V^. 

We wish to solve the equation, 

m-m)m)^iS{l), (3) 

where the underlines indicate matrices. We require the Green's function K(2, 1) to be a 
solution of 

K(2, l)(i7- "^i + m) = (4) 
where the arrow indicates differentiation to the left. 

Our problem is to solve Eq. (3) for ^ at a point 2 given its values at points numbered 
1 on a surface enclosing it and the source S_{1) lying inside the surface. See Figure 1. 
To avoid any singularity at point 2, we introduce a small sphere with radius e, centered 
at the point 2, and shift point 2 to its surface. The radius, e, will be made vanishingly 
small when all other calculations have been finished. To obtain an integral equation, we 
multiply Eq. (3) by K^(2, 1) from the left and integrate xi over the region exterior to the 
sphere containing point 2 and within the outer surface. Then we multiply Eq. (4) from 
the right by ^ and again limit it to the region exterior to the sphere and within the outer 
surface. Finally we add the altered Eq. (4) to the altered Eq. (3) and divide by i. The 
result is 

j K{2, 1)7 • ("^1 + tiMl)d^xi = J K{2, l)S{l)d'^xi (5). 
The left hand side of this equation is equivalent to 

di^ ■ [m, l)7^*(l)]d4a;i = ^[K(2, l)^^'m]d^l^. 
= [ [K{2, l)7'^*(l)]dEi^ - [ [K{2, (6) 

JEq J(T 



where we have made use of Gauss's theorem, and where S is the surface in four dimensions 
that limits the integration. It is made up of the outer surface Eg and the infinitesimal 
spherical surface a containing the point 2. The derivation of the surface integral from the 
volume integral usually gives an outward pointing surface normal 4-vector for dT, and this 
is the case for Eg, but we reverse the inward pointing surface normal 4-vector on a so that 
the surface normal 4-vector da on this infinitesimal spherical surface points outward. This 
has reversed the sign in front of the integral over a in the last equation. 
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When the radius e gets very small, the point 1 in the integral over cr approaches 2 so 
that \1/ can be factored out of it. Replacing the left hand side of Eq. (5) by the right hand 
side of Eq. (6) then gives 



[K(2, l)^^'m)]d^ln - 

So 



J K{2,l)dai^ 7fSL(2) = J Ki2,l)Sil)dW (7) 



In order to determine K_{2, 1), it is convenient to write it in terms of a scalar function 
so that 

K(2,l) = ,^(i?2i)(-i7- d2 + m). (8) 

Then if we multiply this equation from the right by (i'j ■ *d2 + rn), Eq. (4) shows that 
satisfies 

(□i-m2)<^(i?2l) = (9) 

This is the equation satisfied by in the region between the surface a and the surface So. 
The integral over a in Eq. (7) is then, according to Eq (8), 

[ K{2, l)(ic7f = -i7^ [ dir^<j){e)da'( + m0(e) I da^. (10) 
Ja Ja Ja 

The last integral vanishes because da^^ is a component of a vector of constant length that 
will point in every direction with equal weight. Since 

dlr^(f>{R2l) = Xir^4>'{R2i)/R21, (11) 

Eq. (10) reduces to 

[ K{2,l)da^^ = -i-f1 I xi^da'(<p'{e)/e. (12) 
J a J a 

Since u is a four-dimensional sphere, the value of the last integral is the same for 
any value of i]. Let us choose four for i]. The integral over the spatial dimensions will 
then vanish for values of // unequal to four. We shall use 9 and for the usual angles in 
three dimensions and add x cts the angle between an arbitrary vector and x'^. (See Ref. 4, 
Appendix B, for angles in four dimensions.) For a point on cr, 

4 4 
X4 — g^^x — X — e cosx 

and 

AO O 

da^ = e sinO sin x cosx dO d4> dx ■ 
From these two equations, it is clear that 

4d(7t = eV/2 (13) 

'cr 

and 
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xi^da^ = g!^,e^TT'^/2. (14) 
We can substitute this expression into Eq. (12) to obtain 

K(2, 1) daf = -Z7f 6^ (7rV2)0'(6). (15) 
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If we reduce Eq. (9) to a radial equation, it will be Bessel's equation. We need a 
solution with the property that, when substituted into Eq. (12) along with the value of 

the integral in Eq. (14), it will cancel and leave a constant. Such a function is 

(piR) = aH[^\imR)/R = aHf\-imR) / R (16) 

where is a Hankel function. The Bessel equation is linear so that a is an arbitrary- 
constant. This function (j){e) has the limiting form — 2o/(7rme^) as e vanishes. We now 
substitute this into Eq. (15) and the result, in turn, into Eq. (7). Then this equation can 
be rearranged to give 

*(2) = - j K(2, 1)5(1) + ^ K(2, l)7^*(l)dEi^ (17) 

where the constant a has been assigned the value im? /'&tt. 

If we have a boundary So in the form of an initial flat constant time surface, a large 
sphere in space at intermediate times, and a final flat constant time surface, we can undo 
the Wick rotation and obtain the result, 

*(2) = j K{2, l)S{l)d'^xi - y^^(2: lh^-^il)d^xi + J^K{2, l)^''i;{l)d^xi, (18) 

where the subscript i indicates the initial and / indicates the final times. This equation 
agrees with Feynman's"^ Eq.(19) (on his page 754) if there is no source S_. However, 
Feynman's Green's function is given in Eqs. (1) and (2), and ours replaces Eq.(2) by 

/+(2, 1) = 0(i?2l) = {m/8n)H[^\-imR2,i)/R^2,iy (19) 

This differs from Feynman's by only a delta function. A similar derivation can be car- 
ried out for the photon, and, again, the Green's function has no delta function although 
Feynman's does. As Feynman points out, the overlapping of these delta functions produces 
many, perhaps all, of the infinities in QED. We will see in the next section that our Green's 
function for fermions gives the single loop vacuum polarization for the hydrogen lamb shift 
without infinities. 



III. THE VACUUM POLARIZATION FOR THE LAMB SHIFT 



A large number of text books present the calculation of the vacuum polarization con- 
tribution to the Lamb shift, and we need not repeat all of their analysis. The portion that 
interests us is the single- loop contribution to II^'^ that is represented by Fig. 2. It is a loop 
involving an electron and a positron propagator connecting a bare and a dressed vertex. In 
this case, the excluded part of the volume is enclosed by an infinitesimal sphere centered 
on one of the vertices in space-time. Its radius e is made vanishingly small at the end of 
the calculation. The setting up of the expressions to be evaluated would be most direct in 
configuration space-time since the inner excluded region is a sphere in that space. However 
these expressions are, as usual, most easy to evaluate in four-momentum space. For this 
purpose, we can use Feynman's propagator and then subtract out the contributions that 
originate from the excluded sphere. This excludes the delta functions since they lie inside 
the spherical surface. 
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The expression (to lowest order) to be evaluated in this manner is^ 



(20) 



where an overline indicates a four vector, Tr indicates a trace, and S_^^{p) = ^ — m. 
The four-vector q is the momentum of the photon connected to each of the vertices. Our 
first step is to perform the indicated trace. The resulting integral has an integrand with 
two nonidentical factors in its denominators. They can be made identical with the aid of 
Feynman's transforms^. The results can finally be reduced to 



n^-(g) 



2 ^'^ 
m a ' 



where 



and 

In addition, 



27r 7-1 



{g^'^[{l - ry) M{r^) + i^L{r^)] + {q^'q^/q^) 2r^M{r^)] d(3. (21) 



r] = q\l-p'^)/W, 
M{r]) = Mo + Mi{r]), 

L{v) = LQ+r]Li+ L2{r]). 

."^^ u du 

Lo = 
Li = -2 



+ 
°° u^du 
(1+^' 

°° V? du 



Miiv) = 



1 



(1+^^)^ 
1 



(l+?7 + «)2 (l+w)2 



udu = —ln{l + 77), 



(22) 
(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



L2{V) = 



+ 



2r] 



(l + rj + u)'^ (1 + ^)2 (l + u)^ 



du 



= 2[{l + n)lnil + v)-v]- 



(29) 



The integrals Mq, Lq, and Li are independent of ry and, therefore, of q^. When they 
are substituted into 11^^ and multiply another term that is independent of rj, the Fourier 
transform to space-time produces a delta function located at the center of the spherical 
excluded region. For this reason, they should not have been included in the expression for 
n^^^. Similarly, when Mq, Li, or Mq multiply the q^'"' in ry", this produces a D^"^ operating 
on a delta function. Thus when they are transformed back to space-time, they are also 
confined to the center of the sphere. They, too, should not have been included in 11^'^. 
The same is true for the for the term 2r] that appears in L2. This leaves only Mi in Eq(28) 
and the 2(1 -|- rj) ln{l + rf) portion of L2 to be substituted into n^*^. Substituting these 
symbols into then gives 

n^" = ^ /_'^{^''"[-(i-^)Mi+^)+[(i+^)Mi+^)-^]]-^2r7Mi+^)}rf/5- (30) 
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The T] should not be shown because it produces the derivative of a delta function in the 
sphere, and the remainder of the integrand cancels down to leave 

n^"=^ l\g^--'^^\ln{l+r^)d(3. (31) 

Since r] depends upon the integral can be written as twice the integral from zero to 
one. 

In order to make a comparison with the work of others, it is convenient to define x so 
that, from Eq. (22), 

r; =-^(1-/32) = 4-^(1 -x). (32) 
from which (5=1 — 2x. Substituting these expressions into Eq. (31) then gives 

where 

1 2 
T,(q^) = — f x{l-x)ln{l^^x{l-x)]dx. (34) 

To convert the fine structure constant a to where e is the charge associated with each 
vertex, we divide by 47r, with the result that our equation becomes 

'^{'l ) = — 9 / ^(1 — x\ln\\A 9a;(l — x) ) dx. (35) 

27r^ Jo ^ ' 

This equation is identical to Weinberg's^ Eq (11.2.22) obtained by using dimensional renor- 
malization. He uses it to derive a contribution of -27.13MHz to the Lamb shift of hydrogen. 
This number is in agreement with experiment. 

IV. DISCUSSION 



The method presented here for the evaluation of the one loop vacuum polarization 
contribution to the photon self energy avoids any infinities and, therefore, has no need for 
renormalization. The same is true for the self energy diagram for the electron presented 
in Reference 2. It is, therefore possible to compute the Lamb shift for the hydrogen atom 
to the lowest order without renormalization. If such finite techniques can be shown to 
work to all orders, it will indicate that the equations of quantum electrodynamics are 
without defects and that the presence of infinities in calculations with them are due to 
faulty techniques used to solve them. 
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Figure 1 The outer surface Eg includes tj, ty, and tlie space cylindrical surface. The inner 
sphere a has a radius e that goes to zero in the end. The point S lies in the source. 
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Figure 2 Vacuum Polarization Feynman diagram with surface a enclosing the excluded region. 
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